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Abstract 

This work presents a hyper-viscoelastic model, based on the Hencky-logarithmic strain tensor to model 
the response of a Tire Derived Material (TDM) undergoing moderately large deformations. TDM is a 
composite made by cold forging a mix of rubber fibers and grains, obtained by grinding scrap tires, and 
polyurethane binder. The mechanical properties are highly influenced by the presence of voids associated 
with the granular composition and low tensile strength due to the weak connection at the grain-matrix 
interface. For these reasons, TDM use is restricted to applications concerning a limited range of deformations. 
Experimental tests show that a central feature of the response is connected to highly nonlinear behavior of 
the material under volumetric deformation which conventional hyperelastic models fail in predicting. The 
strain energy function presented here is a variant of the exponentiated Hencky strain energy, which for 
moderate strains is as good as the quadratic Hencky model and in the large strain region improves several 
important features from a mathematical point of view. The proposed form of the exponentiated Hencky 
energy possesses a set of parameters uniquely determined in the infinitesimal strain regime and an orthogonal 
set of parameters to determine the nonlinear response. The hyperelastic model is additionally incorporated 
in a finite deformation viscoelasticity framework that accounts for the two main dissipation mechanisms in 
TDMs, one at the microscale level and one at the macroscale level. The new model is capable of predicting 
different deformation modes in a certain range of frequency and amplitude with a unique set of parameters 
with most of them having a clear physical meaning. This translates into an important advantage with 
respect to overcoming the difficulties related to finding a unique set of optimal material parameters as are 
usually encountered fitting polynomial forms of strain energies. Moreover, by comparing the predictions 
from the proposed constitutive model with experimental data we conclude that the new constitutive model 
gives accurate prediction. 
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Figure 1: Tire Derived Material. 


1 Introduction 


In spite of a rapid growth in technologies and development, scrap tire disposal is still an important and unresolved 
environmental engineering issue today. One promising component, from the spectrum of proposed solutions, is 
the recycling of tires into engineering materials. In this paper, the mechanical characterization of a new Tire 
Derived Material (TDM) for structural applications is proposed. This TDM is obtained by grinding scrap tires 
and rubber factory leftovers into grains and fibers, these together with a polyurethane binder are first leveled 
by a roller and then pressed together to form TDM pads (Figure]^. The material, made of Styrene-Butadiene 
Rubber (SBR), the most popular rubber in tire production, has a low cost and easy to implement production 
cycle - alternately Ethylene-Propylene Diene Monomer (EPDM), a rubber used for the production of a wide 
variety of seals can be used. TDMs can be made in various densities with different mechanical properties and 
have been used mainly in railway applications for vibration reduction [25]. The usual composition of TDMs 
results in high compressibility and allows for use only in a moderately large range of deformation. 

Unfortunately the common hyperelastic material models, e.g. Arruda-Boyce [2], Mooney-Rivlin EH, and Og¬ 
den [38] models, fail in describing their behavior in different deformation modes with a unique set of parameters 
[26]; see also Appendix]^ It is also noted that fitting experimental data of elastomeric solids to polynomial-like 
strain energy functions is not an easy task and can lead to oscillating functions with parameters that may not 
have physical meaning [39]. Here, a logarithmic measure is used to describe the mechanical behavior of TDMs. 

Logarithmic strain, typically referred as “true strain”, was first applied to elasticity theory by the geologist 
G.E. Becker [4] [35], who was an instructor of mining and metallurgy at Berkeley from 1875 to 1879. However, 
its introduction is often attributed to P. Ludwik [21] , who defined (one-dimensional) logarithmic strain via the 
integral y in order to measure the extension of a rod of length 1. Today, the logarithmic strain tensor is 
also named after H. Hencky [HITS], who used it in his systematic deduction of an idealized elastic law [29]. 
The Hencky strain measure has many interesting properties, one of the most useful is that it allows for the 
full realization of an uncoupled additive split of volumetric and deviatoric deformations at finite strain. The 
elastic law proposed by Hencky, which is in good agreement with experiments for a wide class of materials for 
moderately large deformations, as Anand demonstrated [T], is induced by the so-called quadratic Hencky strain 
energy: 

Wn{F) := Wn{U) := n\\devslogUf + | [tr(logf/)]2. (1) 

As observed by Hsu, Davies and Royles m, the choice of strain measure can facilitate the transference of the 
nonlinearity in the stress-strain response from the strain-to-stress mapping to the strain measure itself; see also 
Sharda [43]. This in part explains the success of the quadratic Hencky strain energy approach. Notwithstanding, 
the Hencky strain energy also possesses interesting intrinsic mathematical properties. In a series of articles 
[331 El sa El], a family of isotropic volumetric-isochoric decoupled strain energies 


WMF) -WMU) := 


M g/c II dev^ logI/||2 y [tr(logt/)]^ •£ F > 0 

k 2k 

+00 if det F < 0 


( 2 ) 


based on the Hencky-logarithmic strain tensor logU were studied. Here /i > 0 is the infinitesimal shear 
modulus, K, = > 0 is the infinitesimal bulk modulus with A the first Lame constant, /c, k are dimensionless 
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parameters, F = \/ip is the gradient of deformation, U = VF'^F is the right stretch tensor and devnlogt/ = 
log U — ^tr(log U) is the deviatoric part of the logarithmic^ strain tensor log U. This family of exponentiated 
Hencky strain energies improves upon the well-known properties of the original Hencky strain energy. In 
particular, it was recently found that the Hencky energy (not the logarithmic strain itself) exhibits a fundamental 
property: by purely differential geometric reasoning, it was shown [30l [28l |36] (see also muni 13) that 

((detF)!/” .^,80(71)) = dist2,,d,R+.F ((detF)V" .F,f) = ^ [tr(logf/)]2 = 1 (log(det 

= dist2^od,SL(n) ( (det F) Vn ’ = II dev„ log [/||2, (3) 


where distgeod is the canonical left invariant geodesic distance on the Lie group GL^(n) and distgeod,SL(n): 
distgeod,R+-;^ denote the corresponding geodesic distances on the Lie groups SL(n) and M+ -1]^, respectively (see 
[28l[36l). Thus ITh and WeH have the attractive feature that the energies are based directly on a geometrically 
intrinsic distance of the deformation gradient to the group of rigid rotations. 

For small elastic strains, WeH approximates the classical quadratic Hencky strain energy Wh, which is not 
everywhere rank-one convex; moreover in [33] , it is also pointed out that the quadratic Hencky energy has some 
other serious shortcomings. These points being more or less well-known, it is clear that there cannot exist a 
general mathematical well-posedness result for the quadratic Hencky model Wh, although an existence proof 
for small loads based on the implicit function theorem is, of course, possible. The use of © allows for the 
retention of the fundamental geometric property © of the original Hencky strain energy, but at the same time 
alleviates some of its mathematical drawbacks: up to moderate strains, for principal stretches G (0.7,1.4), 
the exponentiated Hencky formulation © is de facto as good as the quadratic Hencky model Wh, and in the 
large strain region it improves several important features from a mathematical point of view. The main feature 
is that the exponentiated Hencky energy Q satisfies the Legendre-Hadamard condition (rank-one convexity) 
in planar elasto-statics [33l|22], i.e. for n = 2. In this case, the energy is even polyconvex, which, together with 
a coercivity estimate, allows for the application of classical theorems for the existence of energy minimizers 

[Ml HI]. 

Despite these advantages, some aspects of the three-dimensional description remain open, since the formu¬ 
lation is not globally rank-one convex. However, in the three-dimensional case, a loss of ellipticity only occurs 
for extreme distortional strains m- This suggests that the exponentiated Hencky energy © retains its full 
suitability for materials that undergo additional (typically irreversible) phenomena based on distortional criteria 
of the Huber-Hencky-VOn-Mises type, as the involved elasticity tensors can thereby be prevented from reaching 
the non-elliptic domain. This is in sharp contrast to the loss of ellipticity of the quadratic Hencky energy Wh, 
which is not related to the distortional energy alone. 

Beside the above mathematical advantages, the exponentiated Hencky energy satisfies a number of additional 
desirable constitutive properties [33]: for example, planar pure Cauchy shear stress always induces biaxial pure 
shear strain; the limit case k -hoc or, equivalently, u = ^ for the linear Poisson’s ratio z/, corresponds to exact 
finite incompressibility; and there exists a certain three parameter subset {k = | /c) such that uniaxial tension 
leads to no lateral contraction if and only if = 0 (i.e. ZD = |/i), as in linear elasticity (see [24] for further 
discussion). Like the quadratic Hencky energy, the exponentiated Hencky energy also satisfies a weakened 
version of Truesdell’s empirical inequalities [7]. 

In this paper, a variation to the volumetric part of © is proposed to capture the high nonlinearity of TDMs 
when subjected to volumetric deformation: 


f^eHm(i") -W,Um{U) : = 


^ g/c||devnlogZ7||^ gfc[tr(logZ7)]2 |tr(logZ7)P •£ ^ > Q, 

k 2k mk 

+OC if det F < 0, 


(4) 


where z^i is the value of the bulk modulus for large deformations and m and k are dimensionless parameters. 
The main advantage of using the modified exponentiated-Hencky energy comes from the fact that the shear 
and bulk modulus are already uniquely determined in the infinitesimal strain regime, while z^i determines the 
nonlinear response, without interfering with p and k. For the modified exponentiated-Hencky strain energy 

^Here and throughout, log denotes the natural logarithm. 
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proposed in this work, the Kirchhoff stress tensor is given by: 

r = Aog[/WeHm(logC/) = 2 //e''■ devg log [/ 


k [tr(log t/)]"ti.(log U) + Kl 

tr(log U) 




while the Cauchy stress tensor is: 

cr = ^)r = 2 /i log t/|p-tr(iog u) . jj 


+ 


^ [tr(log t/)]^-tr(log jj\ |tr(log 17 )P-tr(log U) |tr(log U)\ 

tr(log U) 


• 11 ^. 


( 5 ) 


(6) 


2 Rate-independent response 

In this section we discuss the rate-independent response of the material, which we will refer to as the equilibrium 
response. The physical properties of the TDMs are greatly influenced by the technologies used in manufacturing 
them. Tests have shown that the density and the mixture composition of the material are the parameters 
that most strongly affect its mechanical properties. In total we consider three types of TDMs with the same 
composition but different densities (Table Q and three modes of deformation: shear, uniaxial compression, 
and pseudo-hydrostatic compression. The TDMs studied where made from tires that were shredded into chips, 
mostly 50 mm in size using a rotary shear shredder with two counter-rotating shafts; tire chips were reduced 
to a size smaller than 10 mm in a granulator while most of the steel cords were liberated by a combination 
of shaking screens and wind shifters; the styrene-butadiene rubber (SBR) granules were selected according to 
their dimensions to fit the desired design mix; polyurethane binder was added to the rubber granules mix and 
uniformly distributed; lastly pads of required size and shape were obtained by pressing the compound in a 
mold. In the sections describing the tests, predictions of the model are also shown, all using a fixed set of fitted 
parameters which are discussed in Section 


Tabie 1: Tire Derived Materiai description. 


Material 

Density 
(kg • in-3) 

Composition 

TDM 500 

500 

90% SBR fibers 

TDM 600 

600 

+ 

TDM 800 

800 

10% SBR grains 


2.1 Simple shear 


Shear tests were carried out at Tun Abdul Razak Research Centre (TARRC) in Hertford (UK). The sampies were 
tested with the ciassicai duai iap simpie shear test configuration commoniy used in the tire industry. Sampies 
of 90 mm in width, 50 mm in iength and 20 mm in thickness, were sheared to a shear strain ampiitude of 100% 
of the initiai rubber thickness at the (slow) strain rate of 0.0067 s“^. The procedure is explained in detail in 
[26] . In simple shear the direction of applied displacements does not coincide with the direction of the principal 
stretches; rat her it involves a rotation of axes. The polar decomposition oi F = R'U gives the right Biot stretch 
tensor U = V F of the deformation and the orthogonal polar factor R\ 


U = 


1 

^7^ +4 


■ 2 7 

0 

\ 

2 

7 

0 

7 7 ^ + 2 

0 

, R = ^ - 

-7 

2 

0 

. 0 0 

TTTi J 

+4 

0 

0 

77T4 J 


U can be orthogonally diagonalized to show: 


logU = 


1 

■yy+4 


—7 log Ai 2 log Ai 0 
2 log Ai 7 log Ai 0 
0 0 0 


(7) 
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where ^1 = 2 v 7^ + ^ + 7j eigenvalue of U. 

for this reason detF = 1 and tr(log[/) = 0. The non-zero 
given by: 


Simple shear does not involve a change in volume; 
Kirchhoff stress component ri 2 from equation <§ is 


ri2 = 4/ie 


2 k log 


- (-1/72+4+7) 


log 


^ (^7^ + 4 + 7 ) 


v/7"+4 


(9) 





(a) TDM 500 (b) TDM 600 (c) TDM 800 

Figure 2: Comparison between shear stress corresponding to exponentiated Hencky energy IFeHnu equation 
and experimental tests for different densities. 

Figure shows the ability of the model to capture the shear behavior out to a shear strain of 100%. It is 
to be noted that the TDM 800 sample, Figure |^c), physically failed in the experiment due to cracking and 
crumbling. Thus the poor correlation in Figure T^c) does not reflect poorly on the model. Unfortunately, an 
intact test is not available for TDM 800 in this configuration. Notwithstanding, given the reasonable agreement 
seen in Figures [^a) and[^b), we feel the model performs well in shear. This is in agreement with the findings 
in [33] when WeH was applied to the rubber data of [4^ [18] . 

2.2 Uniaxial compression 

Uniaxial compression tests were performed using a multi-step relaxation procedure. Thin Teflon sheets with 
lubricant were placed between platens and specimen surfaces. Specimen were cylinders with diameter of ~ 27 
mm and length of ~ 15 mm die-cut from a sheet stock. At each step of the loading process, the strain level is 
increased by 5% at a strain rate of e = 0.01 s“^ up to 70% strain. Between each loading step there is a 600 s 
dwell to allow for relaxation of the material (Figure]^. We consider the value of the stress at the end of each 
dwell-interval as the equilibrium stress; these are shown as the red curves in Figs. [^b)-[^d). Note that we take 
the 1-direction to be the axis of compression. 
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(a) Strain history 


(b) TDM 500 




(c) TDM 600 


(d) TDM 800 


Figure 3: Compression test procedure. First figure (a) shows strain history. Blue curve represent the true data. 
Red curve represents the assumed equilibrium response from the data. 


In order to compare the model to the compression data, one requires knowledge of the material’s three- 
dimensional state of deformation. Since the TDMs are compressible we can not make the usual rubber elasticity 
assumptions and require information on the materials transverse response. To evaluate the transverse behavior, 
pictures were taken at the end of every relaxation period with a digital camera mounted on a tripod; see Figure 
The digital images were processed using the image processing toolbox in MATLAB [23]. The white area A, 




35% 



[ 


! 


t i 


70% 


Figure 4: Original image from the digital camera (left column) versus image after processing (right column) for 
TDM 600 at 0%, 35%, 70% defomation. 

measured in pixels in Figure]^ is the region occupied by the TDM sample, and was determined for each image. 
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The height h of the sample is known at each step of the test from the experimental controller. TDMs due to 
their composition are very difficult to cut and usually the samples do not have a straight edge. For this reason, 
we compute an average width of the sample as: 


^=h 


and the average lateral strain as: 


log A 2 = log Ag = log ^ 
Wo 


( 10 ) 


( 11 ) 


where is the lateral dimension at step n and wq is the lateral dimension at the undeformed state. In order 
to appreciate the non-linear compressibility of TDMs, we define, among several possibilities [191 SSl US] , the 
non-linear Poisson’s coefficient z/ as the negative ratio of the lateral extension and axial contraction measured 
in the logarithmic strain: 

( 12 ) 


ly := 


log Ai 


The measured values of i) are shown in Fig. and display the material’s distinct nonlinearity. 

From equation (H- if we consider 5 the value of the uniaxial Cauchy stress, by projection on the Lie-algebra 
J5t(n) of trace-free tensors, we have 


2 , 


2 /i e'' II log[/||"-tr(log U) iQg u ^ ^ ^ 




0 

1 

0 


\ 


0 

1 

" 3 ^ 


(13) 


which leads to the requirement that under uniaxial stress, U has the following form m- 

0 0 \ / e“ 0 0 


U = 


0 


0 

0—0 <3-+ O X 


= 


0 

Q— 2 CL 

0 


0 

o— ia 


This in turn leads to the stress expression 




(14) 


(15) 


where a = | (log Ai — log A 2 ) and x = log Ai + 21 og A 2 are the experimentally known measures of the deformation. 
Using the measured values of a and x we can compare the model’s predicted stress response from equation (15) 
to the measured stresses. Note that one can also project onto the spherical part of equation which gives 


5 = 3 


Ke*’* x + Kie''!*!" 


(16) 


The comparison of equation (15) and (16) to the experimental data is shown in Fig. The results indicate 
good correlation. Note that the material parameter set is the same as used in the comparison of the shear data. 



Figure 5: Non-linear Poisson’s coefficient V evaluated during compression tests. 
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(a) TDM 500 


(b) TDM 600 


(c) TDM 800 


Figure 6: Comparison between compression stress corresponding to modified 
equation (15) and (16), and experimental tests for different densities. 


exponentiated Hencky energy 


2.3 Pseudo-hydrostatic compression 

As a third mode of deformation we consider an experiment that for quasi-incompressible materials provides an 
approximation to a hydrostatic compression test. Here we consider a lubricated cylindrical specimen that is 
inserted into a rigid (steel) cavity of the same radius and then axially compressed. During the test a force was 
applied on top of a steel piston at a volume ratio rate 0.0067s“^. The procedure is explained in details in [26] . 
Considering the axis of compression to be the 1-direction, to good approximation this test follows the kinematic 
path 

/ logAi 0 0 \ 

\ogU= 0 0 0, (17) 

\ 0 0 0 / 

where in the experiment logAi is measured. Likewise in the experiment an = sn +p is measured, where sn 
is the axial deviatoric stress component and p is the pressure. The experiment is designed to test the pressure- 
volume relation. The Jacobian of the deformation, det F = Ai, is readily available from the experiment. However 
the pressure is approximated as p ^ an, which is only valid for p ^ Sn. For the present model, under the 
given deformation state. 


Sn 2/ie^3 

P ^gfc(iogAi)2-iogAi log Ai + Ail^-iogAi 

For TDMs in this experiment, this ratio is not small enough at initial levels of deformation to result in a valid 
pressure-volume experiment. For example, for TDM 500, one must have detF ^ (0.85,1.2) for the ratio to take 
on values of less than 0.05. The plot of the data and the model prediction are shown in Figureand show good 
agreement. Note that the plotted pressure is approximated as an in both the model and the experiment for 
full consistency. It is to be noted that our exponentiated Hencky energy is crucial here. The quadratic Hencky 
energy leads to a pressure-volume relation that is not even invertible for detF > e m- The exponentiated 
form itself alleviates this problem. The modified spherical term which we added to the original exponentiated 
Hencky energy allows for the sharp kink in the pressure-volume relation due to void collapse. As with the prior 
deformation modes, the results shown are produced with the exact same set of parameters. 
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detF detF detF 

(a) TDM 500 (b) TDM 600 (c) TDM 800 

Figure 7: Comparison between pseudo-volumetric response corresponding to equation (§ and experimental 
tests for different densities with p ^ an. 


3 Parameter identification equilibrium response 

For the three states of deformation under consideration, we have utilized a single set of parameters per material 
density. The parameter estimation itself was performed using the non-linear least square (NFS) optimization 
method to minimize the residuals. The main difficulty in such a NFS problem is to find a unique set of optimal 
parameters. Several numerical algorithms have been used in the literature to solve NFS problems mm-, they 
are usually a modification of the Newton method and require an initial guess for the solution. The iterative 
technique furnishes an optimal solution when some stopping criteria are met. In this paper we modified the 
function Isqcurvefit in the optimization Toolbox of MATFAB [24] to fit the different sets of data together [39] . 
The initial guess used was the physical estimate of the moduli obtained by previous experiments on the TDMs 
(Table ||); the remaining parameters were initialized to k = k = k = 0 and m = 2. We imposed positive values 
as lower bounds on the parameters; moreover we respect the mathematical bounds on the parameters [33] . 
Table gives the resulting optimized values which were used for modelling the behavior displayed in the prior 
sections. It is useful to note that the fit values for the shear, low strain bulk, and large strain bulk moduli are 
all sensibly close to the original values derived directly from the tangents to the experimental response. The 
only oddity in values occurs for /^i, the high compression bulk contribution which drops in going from TDM 
600 to TDM 800, whereas an increase would intuitively be expected. The cause for this is as yet unexplained 
but could be related to some sort of local material failure. It should be noted that TDM’s, due to their very 
porous structure, possess structural level failure mechanisms that initiate before the usually expected ones, such 
as filler dewetting umin] etc. as is common in filled elastomeric composites. 
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Table 2: Initial guess for parameter identification procedure from tangents to the experimental data*. 


Material 


K 

Ki 


(MPa) 

(MPa) 

(MPa) 

TDM 500 

0.22 

2.40 

297 

TDM 600 

0.31 

2.70 

315 

TDM 800 

0.63 

4.50 

281 


Parameters obtained from m- 


Table 3: Parameters - exponentiated-Hencky energy function. 


WeHm(C/) 


fc||dev3 logUII 

2 K 

H -^ 

2k 

[tr(log U)] 

2 

+ ^ 
m k 

gfc|tr(logU)r 


Material 

n 

k 

K 

k 

Ki 

k 

m 

(MPa) 

(-) 

(MPa) 

(-) 

(MPa) 

(-) 

(-) 

TDM 500 

0.12 

0.59 

1.40 

0.13 

116 

268 

4 

TDM 600 

0.19 

0.39 

2.80 

0.13 

647 

1989 

6 

TDM 800 

0.50 

0.27 

4.40 

0.13 

404 

1353 
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4 Non-equilibrium response 

In experimental investigations, filler-reinforced rubber like SBR, the main component of TDMs, shows many 
nonlinear effects when subjected to dynamic loads. The main ones being the pronounced dependence of the 
material behaviour on the dynamic strain amplitude together with rate dependent response. The first, often 
termed the Payne-effect [40], can be described as a reversible softening with increasing dynamic strain ampli¬ 
tude. To the authors’ knowledge there are still no well-accepted models that incorporate both the Payne-effect 
and rate dependency. Thus our goal in modeling the non-equilibrium behavior will be limited to the rate de¬ 
pendency of the material at fixed frequencies and amplitudes. To that modest end, a finite strain model of 
viscoelasticity is constructed considering the multiplicative decomposition of the deformation gradient F into 
elastic Fe and inelastic Fi parts as proposed by Sidoroff [44]. Here we assume the existence of two viscous 
mechanisms associated to the material: intermolecular resistance at the microscale level and grain interactions 
at the macroscale level. The first is associated with a Maxwell element including a non-linear spring (A) while 
the second is associated to a Maxwell element in which a linear spring is included (B). The choice of modeling 
the interaction between the rubber particle inside the TDM with a linear law is due to the presence of the 
binder at the grain interface. The binder acts as an internal constraint allowing only normal contact interaction 
between the grains. For this reason, both relative rotation and sliding, which are usually found in granular 
materials [sTiin], are not allowed or are negligible between the grains and fibers. A one dimensional rheological 
schematic is presented in Figure 
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EQ 


Figure 8 : One dimensional rheological model for rate-dependent behavior of TDM. 


If we consider for each viscous mechanism a set of internal variables {k = A^B) that can be viewed as 
the deformation gradient associated with each dashpot, then we can write the total free energy as: 




^NEq(^^) + ^NEQ(^f) 5 


(19) 


where represents the strain energy in the equilibrium spring discussed in the previous part of this paper 

and the strain energy in each Maxwell element associated to the “elastic” left Cauchy deformation tensor 

' [^e ^Iso called the Finger tensor. For most polymer based materials, the volumetric deformation 
is purely elastic and the viscous effects are restricted to the isochoric component of the deformation. Following 
this assumption the strain energy for the Maxwell elements can be written as: 


^ lldevalogb^ll^ 


'NEQ(^e) - 
uB 


W^eq{K ) = mb II dev3 log b, 


B\\2 


( 20 ) 

( 21 ) 


The general theory of viscoelasticity at finite strains used in this work follows the developments of miiii]. Here 
we recall only the essential equations. As a consequence of the Clausius-Duhem inequality, the Kirchhoff stress 
is given as 


^EQ 




'nEQ 5 


( 22 ) 


where t^eq ~ Consistent with the Clausius-Duhem inequality, the evolution of is given by: 


where Cy is the Lie derivative of along the velocity field of the material motion, 

[F^]^F-^^ and is an isotropic fourth order fluidity tensor defined as: 


(23) 


C'i' = 


(24) 


Here is the fourth order symmetric identity tensor, while 77 ^ > 0 represents the deviatoric viscosities. In our 
model 77 ^ = 12 s • N • mm“^ and 77 ^ = 1 s • N • mm“^ for all the different densities and for all the testing modes. 
The model presented implicitly defines the total Kirchhoff stress r. The actual use of the model requires the 


solution of the nonlinear relation (23) which we perform using the predictor-corrector method advocated in m 

and mi. 


4.1 Dynamic shear test 

Dynamic shear tests were performed at TARRC using the dual lap set up with samples of the same dimension 
as used in the static tests. For each sample the displacement was driven up to 33% and 100% of the initial 
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thickness and the tests, for each amplitude, were carried out for 10 cycles. Further, each test was conducted at 
two frequencies 0.1 Hz and 1 Hz. 





Figure 9: Comparison between cyclic shear tests (markers) and the viscoelastic model based on the modified 
exponentiated Hencky energy (solid line), equation ( [l^ , for different frequencies at 100% amplitude. 





Figure 10: Comparison between cyclic shear tests (markers) and the viscoelastic model based on the modified 
exponentiated Hencky energy (solid line), equation ( [l^ , for different amplitudes at 1 Hz. 


Figure shows the dynamic response at 100% strain amplitude for two frequencies. Over this range of 
frequencies the material only weakly depends on the strain rate. Figure considers 100% and 33% strain 
amplitude at a loading frequency of 1 Hz. Here one observes a strong amplitude dependent response. Also 
shown in Figs. and are the predictions from fitting the model to the data. The match is seen to be 
quite acceptable but it should be emphasized that due to the Payne-effect the values of the non-equilibrium 
parameters are frequency and amplitude dependent, as discussed more fully below. 


4.2 Dynamic compression test 

Uniaxial compression tests were carried out on a Bose Electroforce machine in a frequency range of 0.1 Hz to 
25 Hz with the same setup as the static tests. The strain history consists of a static pre-strain of 10% and a 


superimposed sinusoidal excitation varying in amplitude in the range of 1% to 20%. Figures pT] and show 
a few tests (data are shown as markers) from the many performed as they are representative of the overall 
abilities of the model. The model parameters used to generate the solid lines are discussed in the next section. 
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(a) TDM 500 


(b) TDM 600 


(c) TDM 800 


Figure 11: Comparison between cyclic compression tests (markers) and the viscoelastic model based on the 
modified exponentiated Hencky energy (solid line), equation ( p!^ , for different frequencies at 20% ampli¬ 

tude. 


Figureshows the steady-state hysteresis curves at frequencies 0.1 Hz, 1 Hz and 25 Hz at a constant strain 
amplitude of 20%. It indicates that the stress increases with increasing frequencies and the material is stiffer 
at higher frequency. The correlation between model and experiment is seen to be good. 



(a) TDM 500 


(b) TDM 600 


(c) TDM 800 


Figure 12: Comparison between cyclic compression tests (markers) and the viscoelastic model based on the 
modified exponentiated Hencky energy VF^Hm (solid line), equation ( p^ , for different amplitudes at 1 Hz. 


Figure shows the steady-state hysteresis curves with two different strain amplitudes, 10% and 20%, at 
1 Hz frequency. These graphs confirm that the material subjected to smaller compressive strain amplitudes is 
stiffer than material subjected to larger compressive strain amplitudes, similar to what was seen in the shear 
tests and the correlation between data and experiment is good. Other specimens showed a similar pattern even 
though they were taken about different mean strains and frequencies. They are not shown in this paper for 
brevity. As part of the compressive strain campaign, we also evaluated the energy dissipated per hysteresis 
cycle as: 

D=(b(TiidXi. (25) 
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Figure 13: Energy dissipation per hysteresis cycle in compression. 

The results, which are shown in Figure display the presence of two transition regions which validates our 
use of two Maxwell elements in our model for TDMs in this range of amplitude and frequency. This in large 
part contributes to the good ability of the model to capture the hysterisis curves shown in Figs. i-EH 

5 Parameter identification viscous model 


Since the response for large deformations is not a perfect sinusoid, the hysteresis cycles are not elliptical. 
Therefore, the classical definition of storage and loss modulus is inapplicable. In this section we take a look 
at the parameters for each dissipation mechanism and show how they vary with frequency and amplitude. 
The parameters associated with mechanism A, which we term the microscale level, show both amplitude and 
frequency dependence (Figure 14|15 ). The parameter (ia decreases with amplitude then it stays constant both 
in frequency and amplitude. It is mainly amplitude dependent reproducing the Payne-effect well known to be 
present in filler-reinforced rubber. The parameter kA captures the frequency dependency of the material and it 
is constant with amplitude. The single parameter associated with the mechanism B varies with the amplitude 
and stays constant with frequency (Figure [T^. 



(a) TDM 500 (b) TDM 600 (c) TDM 800 

Figure 14: Model parameter /ia in the frequency and amplitude range. 
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(a) TDM 500 (b) TDM 600 (c) TDM 800 

Figure 15; Model parameter kA in the frequency and amplitude range. 
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(c) TDM 800 


Figure 16: Model parameter /i^ in the frequency and amplitude range. 


6 Conclusion 

In this study we presented a hyper-visco-elastic constitutive model for TDMs to characterize the response of this 
class of materials under different deformation modes. The new model is based on an exponentiated Hencky strain 
energy that as shown in previous work improves several mathematical properties with respect to the classical 
quadratic Hencky energy function. There are two main advantages of the newly proposed model. The first one 
is its ability to describe different deformation modes with a unique set of parameters in the equilibrium range. 
The second advantage is that most of the parameters have a physical meaning simplifying the parameter fitting 
procedure. An extensive experimental campaign on TDMs was conducted with both static and dynamic tests. 
It was not the objective of this work to provide a model capable of describing dynamic characteristics of TDMs 
that includes simultaneous frequency and amplitude dependent effects, but rather to explore this new material 
and find the simplest model to characterize it for specific frequencies and amplitude. The predicted results are 
in excellent agreement with the presented data and thus give a viable model for engineering applications of 
TDMs. 
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A Classical model fits 


Montella, Calabrese, and Serino |26] attempted to fit classical hyperelastic models to TDM response but showed 
that the fits in general were poor. With our additional data, we find similar (if not worse) results. Using the 
same fitting method as was done with our new exponentiated Hencky model, one has the following results when 
trying to utilize the Arruda-Boyce model [2] , the Mooney-Rivlin model m, and the Ogden model truncated to 
the third order expansion [38] for compression Eig. 17, shear Eig. 18, and pseudo-compression Eig. 19 




(a) TDM 500 (b) TDM 600 (c) TDM 800 


Eigure 17: Comparison between compression stress corresponding to different hyperelastic models and experi¬ 
mental tests for different densities. 
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(a) TDM 500 (b) TDM 600 (c) TDM 800 

Figure 18: Comparison between shear stress corresponding to different hyperelastic models and experimental 
tests for different densities. 





(c) TDM 800 


Figure 19: Comparison between pseudo-volumetric response for different hyperelastic models and experimental 
tests for different densities. 
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